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Abstract 

By interpreting the fusion matrix as an adjacency matrix we associate a loop model 
to every primary operator of a generic conformal field theory. The weight of these loop 
models is given by the quantum dimension of the corresponding primary operator. Using 
the known results for the 0(n) models we establish a relationship between these models 
and SLEs. The method is applied to WZW, c < 1 minimal conformal field theories and 
other coset models. 
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1 Introduction 

The study of statistical mechanics systems related to loop models is interesting both from 
the physical and the mathematical point of views. Most of the statistical models studied in 
physics, from the Ising and the q-state Potts model to more complex vertex models, can be 
represented in terms of loops. The loop representation of the Ising model is very easy to 
understand: loops correspond to domain walls separating regions of opposite magnetization. 
If the Ising model is defined on the infinite plane or in a box with periodic boundary condi- 
tions, then the boundary of domain walls do not have open ends or branch points. The same 
property is true in the more general Potts model if we define domain walls appropriately. 
*e-mail: Rajabpour@to.infn.it 



A fundamental model, starting point in the construction of more complicated loop systems, 
is the 0(n) model and its connection with exactly solvable models led to the introduction of 
powerful integrable model technologies. A first interesting progress in this direction was the 
discovering of a critical point in the 0(n) model made by Nienhuis [I], however the rigorous 
proof of this result is still missing. 

After the advent of conformal field theory (CFT) as a tool to classify two dimensional 
critical phenomena, the main efforts to understand critical loop models have been based on 
the CFT approach. Many critical loop models were defined and studied also with other 
methods, such as the Yang-Baxter equation [2] and non rigorous coulomb gas methods [3j. 
For a review of the coulomb gas method for more complex loop models see [3]. The critical 
points of these models are described by CFTs, the early versions correspond to well-known 
CFTs such as the c < 1 minimal models, but most of recent proposals links to CFTs with 
extra symmetries like the Wess-Zumino-Witten (WZW) models [5]. For recent progress in 
this direction see [6j [7] . In parallel with the progress on CFTs as a purely algebraic method to 
classify critical points of statistical systems Pasquier in [9] proposed the ADE lattice models, 
a generalization of the restricted solid-on-solid models of Andrews, Baxter and Forester [8], 
as the physical candidates for particular CFTs. Pasquier's work is closely related to the 
Cappelli-Itzykson-Zuber ADE classification of CFTs [10] established by noticing that the 
operator content of many CFTs is constrained by the modular invariance of the partition 
function on torus geometry and it is related to Coxeter diagrams. 

The definition of the corresponding loop models and the study of the boundary changing 
operators for the A series was initiated by Saleur and Bauer in [11]. The A series is related 
to dense critical loop models, while the dilute versions were studied later by Kostov [12J. 

Some progress has been made in the study of loops related to the remaining ADE lattice 
systems, but many of their properties are still unknown and in particular the general fused 
ADE models are not well explored. There are at least a couple of reasons that make the study 
of loop models related to statistical systems interesting. Firstly, from the mathematical point 
of view it may give us good candidates for the Schramm Loewner evolution (SLE), a method 
discovered by Schramm [13] to classify conformally invariant curves connecting two distinct 
boundary points in a simply connected domain. The parameter describing the curve is the 
drift parameter k (For recent reviews see [HIE].). The exact relation between SLEs and 
CFTs was discovered by Bauer and Bernard [16]: using the method of null vectors of CFT 
they were able to find a simple relation between the CFT central charge and the drift At. This 
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method, based on the original results of Schramm, was generalized later to loops without 
open ends by Sheffield and Werner [17] and named Conformal Loop Ensemble (CLE). The 
connection between ADE loops and CLE was studied by Cardy [18] in both dilute and dense 
cases by mapping the height variable of ADE lattice models to the 0{n) model. 

Another reason of interest is related to the study of topological quantum field theory in 
2+1 dimensions and to its application to topological quantum computation [19J. The ground- 
state wave function of these topological models coincides with the loop ensemble. This means 
that the ground-state correlators in 2 + 1 dimensional topological quantum field theories are 
equal to particular correlators in 2D classical loop models |20j . The corresponding loops 
have weights related to the quantum dimension of the operators appearing in the model. It 
seems that mathematically one can classify the ground-state of topological theories by two 
dimensional critical loop models or equivalently by CFTs. 

The above considerations suggest that a classification of the loop representations related 
to CFTs may be relevant for a wide range of applications. 

In the following, using the methods introduced by Cardy in [18j . we shall define a loop 
representation for a generic CFT. To define these loop models we use the fusion matrix of 
CFTs as an adjacency matrix and map the resulting height model on a 0(n) system. We 
find a different 0{n) model for every CFT primary operator and observe that the weights 
of the loops coincide with the quantum dimensions of the corresponding primary operator. 
Using the link between 0(n) and CLEs one can map every CFT to SLEs and to CLEs. 

This allows the investigation of the link between SLEs and CFTs to be carried out without 
the use of, the sometimes quite complicated, null vector method. The main ingredient of the 
method is the S matrix or the complete form of fusion matrix which is complex for some 
extended CFTs. 

The paper is organized as follows: in section 2 we use the fusion matrix as an adjacency 
matrix and define a height model using the S matrix as a weight for the plaquettes in a 
triangular lattice. 

Following [18] , using the height model we will associate to every height configuration an 
ensemble of loops, and by the Verlinde formula find that the partition function of the loop 
ensemble matches that of a particular O(n) model. 

In section 3 we investigate the loop representation of some simple CFTs such as the Ising 
and the tricritical Ising models, the c < 1 minimal models, the WZW and in particular the 
SUk(2) models. Our prediction for the SLE drift parameter k is in agreement with all the 
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previously known examples. 

In section 4 we will generalize the methods introduced in section 2 and define loop models 
on the square lattice related to the dense phase of the 0(n) models. The definition of loops on 
the square lattice is more natural because the models arising from our definition are always 
critical without the need of a fine-tuning of the parameters. In section 5 we will briefly discuss 
the connection between these results and boundary CFT and show that the models are well 
defined also in presence of a boundary. 

Section 6 contains our conclusions with a brief description of the work in progress moti- 
vated by these results. 

2 Loop models for general CFTs 



To define loop models for a given CFT we shall use the modular S matrix. To define 
the S matrix we first need the expression of the modular invariant partition function on the 
torus geometry in terms of the characters and the parameter r = e r , where R is the space 
length of the system and (5 is the time period, as follows 

Z = J2 M h,hXh(r)Xh(r)- (2-1) 

In (|2.ip the sum is over the operator content of the theory and M h ^ is the non-negative integer 
that specifies how many times a given representation enters, that is the number of primary 
operators with scaling dimensions (h,h). The function Xh = tr(T L °~2i) is the character of 
the theory with the sum over the descendants of the primary field with conformal weight h. 

From now on we shall use a, b, c to label the primary operators (they are not the weights 
of primary operators) in the theory independent of the holomorphicity of the operator. Then 
one can define the modular S matrix by the modular transformation which does not change 
the partition function but acts nontrivially on the characters 

Xa(~) =E^X&(t). (2.2) 
v ' b 

Another important ingredient is the fusion rule related to the operator product expansion of 
a pair of primary operators, 

4> a 4b = Y. NC oh4>c (2.3) 

c 
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Figure 1: A triangular plaquette with a ^ b and the corresponding curve segment on the 
dual honeycomb lattice. 

In (|2.3p N£ b is the fusion coefficient and the RHS involves the entire tower of descendants 

of the primary operator <p c . It is worthwhile to mention that it is possible to interpret the 

fusion coefficients as the elements of a matrix (N a ) b = N£ b , so we have a fusion matrix for 

every primary operator. 

To define a loop model we interpret the fusion matrix as the adjacency matrix associated 

to a graph [18]. The graph of a primary operator <p a has g vertices where g is the number of 

primary operators in the theory and edges connecting pairs of vertices (6, c) when N£ b = 1. 

Following [18J one can define a height model on the triangular lattice by imposing that 

the height hj at the site j can take values 0, 1, . . . ,g — 1. Then constraint the heights at 

neighboring sites according to the incidence matrix associated to a given primary field <p a : 

only neighbor heights hk and hk with {N a )^ = 1 are admissible. One should notice that 

according to the fusion matrix rule neighbor heights maybe identical, which is not the case 

for ADE models. Actually, for a consistent definition of loop models on a triangular lattice 

at least two of the heights at the corners of an elementary triangular plaquette should be 

equal then the weights for the elementary plaquette are defined as follows. If the heights 

of plaquette are (c,b,b) with c ^ b then weight is x(4l) 1//6 , where S is the modular matrix 

and k is arbitrary. If the heights are all equal then the weight is 1 except for those with 

N% b ^ that have weights 1 or x depending on the particular model considered, as it will be 

explained below in the paper. 

The next step is to mark triangles with unequal heights (c, b, b) drawing a curved segment 

on the dual honeycomb lattice [18] and linking to the center the midpoints of the two edges 

1 We can have loops in our graph but at least for minimal models and restricting the analysis to primary 
operators of interest, it is possible to prove that there are not graphs with loops smaller or equal to three. In 
general we just define models for graphs with loops longer than three or even without loops. 
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with different heights (b and c) at the extremes (See figure 1). The difference with the more 
standard ADE models is related to the equal height (b,b,b) plaquette with N% b = 1. For 
plaquettes with height (b, b, b) there are three possibilities for drawing curve segments. In 
these cases we will suppose, as in percolation problems, that the lines in dual honeycomb 
lattice choose randomly two of the edges of triangular lattice consistently with the other sites. 
This define a loop configuration for every height configuration. Summing over the admissible 
values of heights consistent with a given loop configuration we find 

£(«f = f , (2.4) 

where the sum is just over b. To get this formula we have used the Verlinde formula for 

CFTs, X^(iV a )S;T& = which means that the 6th element of the eigenvector of N a 

s a 

with eigenvalue is given by We take always k = to get the largest eigenvalue and 
have positive real weights in our height models. For properly-defined fusion matrices these 
eigenvalues will coincide with the quantum dimension of (p a : 

da = |j. (2.5) 

Summing iteratively over all heights of all clusters gives a factor d a for each closed loop, so 
the partition function of our model has a 0{ri) like partition function 

Z = 5>C (2.6) 

where I is the number of bonds in the loop configuration and N is the number of loops. 

Using this method we are then able to associate to every primary field of a given CFT 
an 0(n) model with n = d a . Let us discuss some simple but general properties of these 
models. Firstly it is not difficult to see that for the identity operator our incidence graph 
has g disconnected points with a blob attached to them to indicate that they are adjacent 
to themselves. This means that the configuration space is decomposable and that all the 
heights should be equal. In loop language means that we have a percolation or an Ising loop 
model depending on the value of x. It is not difficult to see that for any non-connected graph 
our height configuration space is decomposable and that we can treat the different parts 
of our graph separately. Another interesting fact is related to the consistency of different 
parts of graph: if we write for different connected components of graph the corresponding 
adjacency matrices then all the parts will have the same loop weights and the same quantum 
dimension. Therefore, we can consistently keep only one part of the graph and find the loop 



model associated to it. We shall clarify further this issue later by discussing some minimal 
model examples in details. 

It was found long time ago in pQ that the 0(n) model posses a dilute critical point for n < 2 
with x c = -, 1 = : correspondingly our loop models will have a critical point just for the 

y 2+-\/2— n 

fields with quantum dimension smaller than 2. The 0(n) model has another critical regime, 
the so-called dense phase, for x = (x c , oo) which corresponds to a different universality class. 
The mapping to the 0{n) model helps us to find the connection with CLE. From coulomb 
gas arguments we know that, in the dilute regime, the loop weight has the following relation 
with the drift in the CLE equation, 

4-7T 

d a = -2cos(— ). (2.7) 

We will use this equation to establish a connection between our models and CLEs and find 
that the prediction is in excellent agreement with previous results. We should mention that 
for the dense phase the above equation is still true if we work in the region 4 < k < 8. In the 
next section we will investigate the loop representation of some simple CFTs. 



3 Some examples 

In this section we show how the method described in the previous section is working for 
some simple models such as minimal models, SUk(N) WZW models and coset models. 

3.1 Minimal Models: 

Consider the minimal models M(p',p), the simplest case is the Ising model M(4, 3). This 
model has three primary operators [5j: 1,<t, e with graphs depicted in figure 2. We discussed 
in the previous section how the operator 1 is related to percolation or Ising model with k = 6 
or 3. The other two cases are more interesting, for example the graph of N e has the part 
similar to the A2 graph and a disconnected point, using d e = — 2 cos(7r|) = 1, it is not difficult 
to see that both parts have the same dominant eigenvalues corresponding to k = 3 which is 
in agreement with [18]. The graph for N a is similar to A3 and gives us k = ^ in the dense 
phase, which is the dual of the previous case. It is quite surprising that the critical behavior 
of these loops are completely in agreement with the known Ising model. 

The next simple minimal model is tricritical Ising model M(5, 4) with six primary op- 
erators and fifteen different fusion relations [5J. The fusion matrix graphs for this model 
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Figure 2: Graphs of fusion matrices of primary operators in Ising model, from left to right 
the fusion graph of 1, e and a. 

are shown in the Fig 3. The familiar cases with the critical loop models are l,a',e. The 



again in agreement with [18] . The case N e is similar to and related to k = 5. One can 
check simply that the other part of the graph has the same eigenvalue and so related to the 
same kind of loop model, this kind of graphs are the first kind of non-familiar height models 
corresponding to loop models. We also notice that one part of this graph is similar to the 
second level of the fusion of model. We have two more graphs related to critical loop 
models, N e > has quantum dimension equal to the N e so they have the same properties. The 
quantum dimension related to N e n is equal to one and so again we have percolation or Ising 
model. The last one is the N a which has quantum dimension bigger than two and so we can 
not associate a critical loop model to this primary operator, we notice that this graph is also 
similar to ground sate adjacency diagram of the second level fusion of Aq model [6]. 

Let's summarize the meaning of disconnectedness in the graphs of different fusion matrices 
of primary operators. For both of the above models if we treat the different parts of graph 
as the different models then we can find finally the same loop models. This means that 
the different blocks of our block diagonalizable adjacency matrices have the same largest 
eigenvalues and here our definition of loop model is consistent. One can always choose one 
block and get the corresponding loop model. This is exactly the same as the definition of the 
well known ADE models, for these models always getting one block of our fusion matrix is 
enough. The degeneracy in the eigenvalues of the fusion matrix is the natural result comes 
from the most general property of fusion matrix: different fusion matrices commute with 
themselves. It seems that the above argument is a general property of many different CFTs 
specially the more general cases which we will investigate in the following. 



graph of N a t has two ^3 parts, using d a > 
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Figure 3: Graphs of fusion matrices of primary operators in three critical Ising model in the 
upper row from left to right the fusion graph of 1, e and e', in the lower row from left to right 
the fusion graph of e", a and a'. 



The above calculation is tractable in the general case of minimal models M(p',p) with 
Kac formula for the conformal weights of primary operators 

_ (p'r -ps) 2 - (p-p') 2 

where (rs) is the label of primary operators in the Kac table and 1 < r < p— 1, 1 < s < p' — 1. 
For the known primary operators in the non-extended Kac table, it is not difficult to see that 
because of (r, s) = (p — r,p' — s) there is a huge redundancy in this rectangle. The S matrix 
for minimal models has the following form 



Srs, P a = 2 x ^-{-\) 1 + s P+™ sm(Ap) S m(ir?-sa). (3.2) 

y pp p p 

Using equation (|2.5p the quantum dimension of primary operator <j) rs has the following form 

sm(ir^r) sm(ir-s) 

drs = (-D r+s : % . (3.3) 

sm(7r^)sm(7r^) 

We should mention that the above equation respects the symmetry (r, s) = (p — r, p' — s) 
for the unitary minimal models. One can investigate the properties of graphs for different 
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operators but the most interesting cases are (r, s) = (1,2) and (r, s) = (2,1). The fusion 
matrices for these cases are the following symmetric matrices, 

N(™Z)( r ,s) = $rn,r(Sn,s+l + 5n,s-l) ^(™l)(r,s) = ^n,s(^m,r+l + $m,r-l)- (3.4) 

It is not difficult to see that the graph for these two cases do not have any loop and so we 
do not need to worry about possibility of definition of loop model. For example one can 
show that the graph of (r, s) = (1,2) for unitary minimal models, M(p + l,p), has always 
a part similar to the graph A p ^i and (r, s) = (2, 1) has a part similar to A p . The quantum 
dimension of (r, s) = (1, 2) is 

P 

du = — 2cos(-7r — ). (3.5) 
P 

If we suppose that the loop model is at the critical regime then it is easy to see that this loop 
is related to CLE with the following drift, 

* = 4^, (3.6) 
p 

which is exactly the same as the k predicted for minimal models by using null vectors |16j . 
One can see this by rewriting the level 2 null vector of minimal models as follows, 

L ~ 2 ~ 2{2h + l) Ll1 = L - 2 " jK 1 ' (3 - 7) 

for the unitary Virasoro minimal models and comparing it with the null vector like relation in 
SLE [32]. It seems that there is a close relation between this loop model for (r, s) = (1, 2) and 
CLE. The same calculation is possible for (r, s) = (2,1) and the result is c?2i = — 2cos(-7r^) 
with k = 4^- which is apparently the dual of the previous case and consistent with CLE 
prediction^. In this construction it is also possible to get other CLE drifts by considering 
the dense (dilute) phase of the loop model corresponding to 4>\,2 (^2,1 )• I n this level it is not 
clear how one can choose the right phase but the algorithm for minimal models as we noticed 
is taking dilute phase for ^ 2 an d dense phase for 02. i- For exact correspondence one need 
to investigate the problem in the partition function level. 

These two cases are the familiar cases but as we noticed in tricritical Ising model these are 
not the only loop models with critical point, in fact it is possible to have many critical loop 
models for a minimal model with generic p and p' just we need to have quantum dimension 



2 The outer boundary of fractal curve with k > 4 is given by an other fractal curve with k := 16/k named 
the dual curve. This curve has the same central charge as the original one. 
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smaller or equal to two. Appearing the A p -\ and Ap*-i graphs in our fusion matrix is not 
just an accident because in fact these are the graphs that already appeared in the ADE 
classification of minimal models [10J. From the other point we know that these graphs 
have lattice statistical physics counterparts with the loop representations investigated in [18] , 
so appearing a close connection between classification of modular invariant CFTs in two 
dimensions with respect to the ADE graphs with our loop representation is not so much 
surprising. At least for the minimal models this connection is obvious and the fusion graphs 
of some special operators are exactly the same as the ADE diagram of the corresponding 
CFT. 

For later purposes let's write the results explicitly for the case M(6, 5) which is describing 
the three state Potts model, using equation (|3.3p we find the following quantum dimensions, 

7T 7T 

= 1, d(i, 2 ) = 2cos(-), d (lj3) = 2cos(-), d( M ) = 1, 

d (2) i) = VS, d (3)1) =2, d (2j2 ) = 2v / 3cos(^), d (2j3) = 2V3cos(^), 

d{3,3) = 4cos(^), d( 3)2 ) = 4cos(^), d( 2)4 ) = V3, d (3j4) = -2. (3 

It is evident that different primary fields can have the same quantum dimension and so the 
same loop representations. For example «i )2 = «i i3 = -§■ and «2,i = k 2 ,4 = however 
we know from the connection between CFT and SLE that the fields (1,2) and (2,1) are 
connected to SLE by null vector so we need to answer this question how we can choose 
the right operator. For three state Potts model we can answer the question by looking 
at the modular invariancy of the model which is responsible for constraining the operator 
content of the model to just six operators, (1,1), (2,1), (3,1), (4,1), (3,3), (1,3) in the 
spin representation of the model [5]. However for Fortuin-Kasteleyn representation of the 
model this operator content is not enough and we should insert disorder operators as well. 
The theory has extra W symmetry which is responsible in decreasing the primary operators 
of the theory from ten to six. If we look at the model with the above symmetry then we 
can show that this model is not just a sub theory of minimal model M(6,5), the fields in 
this model represent different primary fields with different characters however the conformal 
dimensions are equal. Using operators W(.z) and W(z), with weight equal to three, one 



can show that 0( 1)3 ) = W_iW_i0( 1)2 ), 0(2,3) = W-i/20(2,2) and 0( 4)1 ) = W_i/20(2,i)> where 
W n and W n ; n = 0, ±1, ±2, ... are the mode operators. The above relations show that the 
operators with the same quantum dimension are in fact related with the W symmetry of 
the model. In the other word they are appearing in the modular partition function together 
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[5], so the equality of quantum dimension of different fields are a sign for the internal extra 
symmetry of the model. This interpretation of three states Potts model is also consistent 
with the parafermionic interpretation [21J. The same story is true for tricritical Ising model 
which has also supersymmetry at the critical point. This symmetry is responsible for the 
equality of quantum dimension of the operators I and e" and also the operators e and e'. 
Both couples come from Neveu-Schwarz sector of superconformal symmetry of the model [5] 
and one is the descendent of the other in the presence of supersymmetry. The other important 
thing to mention is to get the connection to the physical system we need to investigate the 
symmetries of the statistical model specially those that are connected to the symmetries of 
the domain walls. For example for three states Potts model the fluctuating domain walls 
should respect the Z$ symmetry. We know from the boundary conformal operators that the 
corresponding responsible operator is the spin operator [22], but free boundary condition is 
not respecting Z% symmetry and so in the modular invariant partition function we do not 
have the corresponding operator, so we always need to have some informations coming from 
the boundary CFT as well. 

The above calculation can be generalized for the extended Kac table too, for example take 
the simple case M(3, 2) which is related to CFT with zero central charge. The loop model 
for the operator (1,2) has k = I which is self avoiding random walk with zero quantum 
dimension. The loop model of (2, 1) has k = 6 which is percolation. These two cases are 
again consistent with the known results. We should mention that the above naive argument 
is just coming from equation (|3.6[) and does not mean that the adjacency matrix is just A\ 
. The other simple case is M{2, 1) with quantum dimension —2. In this case we have n = 2 
and k = 8, corresponding to loop erased random walk and spanning trees respectively. One 
can simply investigate the extended Kac table of more complex CFTs like Ising model and 
find some new loop models with zero or negative quantum dimension which is related to 
0(n) models. For example in Ising model (1,4) has zero quantum dimension and (1,5) has 
negative quantum dimension. The zero and negative quantum dimensions can not be related 
to unitary CFTs and as it is quite well known the extended Kac table models are not unitary, 
they are related to logarithmic CFTs (LCFT) with special kinds of fusion matrices. For 
example let's briefly see the most familiar LCFT, c = —2. For this case if we just suppose 
the fields with (r, s) = (1, 1), (1,2), (1,3), (1,4), (1,5) then the fusions of the operators make 
a closed algebra [23]. Using the fusion matrices of operator (1,2), which has also integer 
elements more than one, will give us the same dominant eigenvalue equal to zero, this is the 
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same as our above naive argument. We think that these loop models are in close relation with 
the lattice logarithmic minimal models introduced by Pearce, Rasmussen and Zuber [24J. 



3.2 SU k (N): 

Let's investigate the CFTs with additional symmetries specially WZW models and for sim- 
plicity the most simple and familiar case i.e, SUj t {2). The fusion matrix of this models are 
well known; for example see [5], by generating loop model as before we will have some loops 
with weights equal to the quantum dimension of the corresponding primary operators, for 
this case quantum dimensions are completely familiar and has the following form 



_ sin(vr(fl + l)/(k + 2)) 
da ~ sin(7r/(A; + 2)) ' ^ 

where a = 2j with j = 0, 1/2, k/2 is the spin of the representation. The only nontrivial 
case with quantum dimension smaller than 2 is related to a = 1, the spin 1/2, which has 
quantum dimension d\ = 2cos(-i^—). The CLE drift is 

. = 4— (3.10) 

the same as the |25j result which was found by using null vector relations. We notice that 
just for k = 4 we have another quantum dimension which is equal to 2 and related to a = 2, 
spin 1, the corresponding CLE drift is k = 4. We notice that this result is the same as the 
result coming from A^+i with the same adjacency graph and so in the close connection with 
the modular invariant partition functions of this model. The other example is SUk(3) with 
the quantum dimensions, 

_ sin(7r(q + l)/(fc + 3)) sin(7r(6 + l)/(fc + 3)) sin(7r(a + b + 2)/{k + 3)) 
a ' 6 ~ sin(7r/(fc + 3))sin(7r/(fe + 3))sin(27r/(A : + 3)) ' 1 ' ' 

where a, b > and a + b < k. For the simplest case k = 1 we have two quantum dimensions 
di t o = 1 and cfo.i = 1, for k = 2 we have three different quantum dimensions di t o = do,i = 
d\ t \ = 2 cos(|), ^2,0 = ^0,2 = 1- It seems that there is not a nice compact form for the general 
case but the calculation for larger ks is straightforward and shows that most of the quantum 
dimensions are greater than two, except some operators with quantum dimensions equal to 
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one. The same calculation for the SUk(N) is tractable and for k = 2 one of the critical loop 
models has the following CLE drift, 

« (3-12) 

This is similar to equation (|3.1(jp if we replace k with N, this is a general property in SUk(N) 
models and reminiscent of level-rank duality. One can find the quantum dimension of the 
most general affine Lie algebras as the weights of loop models and try to find the CLE drift 
by using 0(n) results, to find quantum dimensions of many CFT's see [5]. 

3.3 Coset Models 



Z{k) model: The same calculation can be done for the Z{k) models related to the coset, 



SU k (2) 



, with the same fusion rules and operator content but with the different conformal 



weights. The central charge of this model is c = 2 an d the conformal weights of the 
primary fields 4>i^ m are given by A^ m = — ijfc with < / < k, < \m\ < 2k — 1 

and I — m £ 2Z so that we should just take half of the grid to remove the redundancy, the 
spin operator is corresponding to m = I in this notation. The modular S matrix are given 
explicitly by 

Sf 1 - = ; ; sin v , JKJ ' exp(i-7-.)- ( 3 - 13 ) 

l ' 3 y/k{k + 2) k + 2 yy k 1 y J 

Since these models have the same fusion rules as the SUk{2) models we find the same quantum 
dimensions and so re-derive the equation (|3.10j) for the first spin operator of this model. For 
this case if we go to the dense phase of the loop model then we have 

« = 4£±f, (3.14) 

which is the same as the formula proposed in [26J for the lattice Z{k) models in the FK rep- 
resentation. We were not able so far to find the formula for the spin representation of these 
models by the familiar S matrix of Z{k) models. The Z{2) case is the Ising model and it is 
equal to the minimal model M(4, 3) and so there is not any ambiguity in this case. For Z(3) 
it is true that the model is equal to three states Potts model in the spin representation how- 
ever the operator content of this model is not equal to M(6, 5). However the two disordered 
operators appearing in the Dihedral description of the 3 state Potts model is missing in this 
description but the equation (|3,10p is giving true answer for k = 3. For higher fe's the result 
is different from the result of [26J, we will come back again to this problem when we shall 
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discuss loop models on square lattice. The above argument means that we should be careful 
in extending our results to the loop models in the lattice Z(k) model defined by spin variables. 



SUk SU^(2) ^ : ^ ne similarity of the quantum dimension of S tfffi and SUk(2) comes from 
the general property of coset models, however coset models can have the larger operator 
content with different conformal weights but they have most of the quantum dimensions of 
original affine algebra. For example take the coset, g ^g^+K2) ^ ' ^ or ^ = 1 h is equal to 
M(p + l,p) minimal models with k = 4|±| for operator (1,2) which is equal to the SUk(2) 
case however the M(p + l,p) minimal models have other critical loop representations as well. 
Generalization of these results to more complicated coset models is straightforward and we 
discussed some of them in the Appendix. 



4 Loop Models on Square Lattice 



The definition of loop models on square lattice is similar to the honeycomb lattice. From 
some point of views it is more natural, for example for the honeycomb lattice we introduced a 
non-known parameter x which it should be fine tuned to get a critical loop model but on the 
square lattice just the S matrix is enough to get a critical loop model. The definition is as 
follows: put some heights on the vertices of the lattice so that the neighbor sites be adjacent 
on the fusion matrix graph. As the honeycomb lattice if the fusion graph has no cycles 
smaller than five then with probability one, one of the diagonally opposite pairs of heights in 
each elementary plaquette should be equal, we will use this constraint for definition of loops. 
Labeling the heights by (a,b,c,d), the weight around elementary plaquette is 



W(a,b,c,d) = (^f )*8 ac + \(^r)*S bd , (4.1) 

where A = 1 for the unitary minimal models. So far the definition is similar to the non- 
dilute ADE models but as we noticed in the first section there are some fusion matrices with 
diagonal terms. This means that it is possible to have some neighboring vertices with the 
same heights, but still we will have at least a diagonally opposite pairs of heights in each 
elementary plaquette with the same height for most of the cases that we are interested. For 
these cases we define the weight of the plaquette by using loop model as follows: in each 
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Figure 4: A plaquette of square lattice with a = c and the corresponding curve segment on 
the dual square lattice 

elementary plaquette we draw an edge connecting the diagonal sites if they are equal then 
we will have a Fortuin-Kasteleyn like graph [18] and then the definition of loop model on the 
medial lattice is straightforward just as Pott's model, Fig 4. Each part of the loop inside 
plaquette carries a weight, (^-j . For plaquette with three equal heights for drawing the 
edges we do not have any ambiguity and we can use the equation (14, ip to get the weight of 
the plaquette. For plaquette with four equal heights we chose one of the edges like percolation 
with weights equal to one and after that the definition of loop model is straightforward. If we 
consider the weights of all parts of loop and then sum over all values of the heights consistent 
with a configuration of loops we will have an 0(n) model with x = 1 which is exactly the 
same as the loops in the FK representation of the Q = d\ state Potts model. Following |31] 
there should be some holomorphic operators with conformal weights equal to s satisfying the 
equation d a = 2sin(s^). This holomorphic operators describe the probability that the curve 
passes between two close points. Then the corresponding drift term of CLE is 



The above argument shows that we can map all of the loop models defined by the above 
method to the Potts model and then by using the known results of critical Potts model we 
can get all of the properties of loop models. Let's investigate some examples: 

Example 1: Minimal models 
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For the case (r, s) = (2, 1) it is easy to see that the spin of corresponding operator should 
be equal to hs^i = 2y — 1 and we have k = 4^- consistent with the argument in [31]. One can 
find the results for the other cases by using the same technics. In all of the cases the results 
are consistent with the results of the previous section. 



Example 2: Z(k) model 

The calculation for the parafermionic models are the same and we can find easily the 
connection of these CFTs to CLE by the above method. One of the famous examples is Z(k) 
model with the following connection to CLE for the operators with spin ^ representation, 



«-4*±f (4.3) 

where the corresponding holomorphic operator describing the probability that the curve 
passes between two close points is, s = We conjecture that the above formula, as 

Santachiara pointed out, could correspond to the FK representation of lattice Z(k) model 
with k > 3. However the result is coming again from the different operator, in the author 
found the above equation by using one of the disordered operators comes from dihedral group 
but we do not have these operators in our calculations. This is another example for different 
primary operators with the same quantum dimensions. We will face again with this ambi- 
guity when we discuss lattice Z(4) model, Ashkin- Teller, as an orbifold model [32J. However 
there are some evidences that Z{k) parafermionic model is a good candidate for the lattice 
Z(k) model [33] but it seems that it is not enough for describing all of the properties of the 
model. In addition it shows that it is important to investigate the fusion matrix of the dihe- 
dral group to improve our understanding about lattice Z(k) models. For Z(k) parafermionic 
CFT there is another more straightforward evidence that the formula (|4.4p is a good candi- 
date, the argument is as follows: In equation (|4.1|) if we put A = — 1 then the corresponding 
CFT describing the height model is Z(k) model which is also antiferromagnetic Potts model 
[7J; see also |34| . and related to loop model with weight — d a . then again by following [31] 
it is easy to show that there is antiholomorphic operator with spin d a = — 2sin(s : |), this 
antiholomorphic operators are describing the probability that the curve passes between two 
close points but in this time by the following connection to CLE drift 

(4.4) 
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which again gives the equation (|4.4p . The tricky point is we need to work with the antiholo- 
morphic operator in this case. 

To close this section we just should mention that the loop models coming from the defi- 
nition on the square lattice is just the dense phase of loop models discussed in the previous 
sections, this is easy to understand if we suppose x = 1 in the 0(n) model's partition function. 
This is always true because for — 2 < n < 2, x c is always smaller than one. 



5 Discussion: connection to BCFT 



In this section we would like to discuss briefly the connection of the argument of section 2 
to boundary conformal field theory (BCFT) and SLE. Firstly let's define boundary condition 
for height models by following 118]: For the height models on the square lattice suppose 
the wired boundary conditions on the real line but with height a on the left part and b on 
the right side. If we work on the half plane then as well as some nested loops we will have a 
curve going from the origin to infinity if a and b are adjacent on the corresponding graph, if 
they don't then it is possible to have more than one curve. Let's first suppose the adjacent 
case then we will have a curve and following [31] it is possible to relate the probability that 
the curve passes between two closed points to a holomorphic operator with spin s, mentioned 
in the previous section, with the expectation value < ip(z) >~ -p. This argument could be 
the starting point for a rigorous proof of the connection of these curves to SLE as Smirnov 
used it to prove the convergence in the Ising case |35j . The above definition is possible for 
all of the examples that we discussed in section 3 independent of the null vector property of 
the corresponding primary operator. This is paradoxical because we know from the work of 
Bauer and Bernard [16], and Friedrich and Werner [36] that the existence of null vectors at 
level two is essential to get SLE. The generalized stochastic equations of [25, 26J showed that 
the extra symmetries in the theory can modify the level two null vector to more complicated 
null vector, so the Virasoro null vector is not complete for describing theories with extra 
symmetries. This can be one of the reasons that we can get a large set of loop models for 
different CFTs without concerning the null vector properties of operators, however the exact 
connection is missing so far. 

The other important point in the connection of our loop models to BCFT is: why we are 
just working with bulk operators without speaking about boundary states? This comes from 
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the old result of Cardy [37] stating that there is a bijection between the possible conformally 
invariant boundary conditions and the bulk primary operators. Let's see this more carefully: 
in the upper half plane for the above boundary condition it was shown in [37J that only one 
copy of the chiral algebra acts and so it is possible to write the partition function on the 
cylinder as a sum on the characters of chiral parts 



where sum is over all of the primary operators and n? b is the element of matrix n.j satisfying 
Verlinde fusion algebra [38] with the eigenvalues the same as for the fusion matrix. So one can 
write the graph corresponding to rtj and find the same loop representation as we found. This 
is the method followed by [38] to classify all the conformal boundary conditions of rational 
CFTs specially SU(2) WZW models. The above argument teaches us that the discussed 
loop models are compatible with conformal boundary conditions and so with SLE however it 
seems that this statement is highly nontrivial and need more investigation [39|. For example 
however we used the Cardy's equation as a consistency condition for conformal boundary 
conditions but it is not trivial why in the lattice model when we go to the continuum the 
model is conformally invariant. This comes from an old problem in statistical physics at the 
critical point that: is the statistical mechanics models at this point conformally invariant 
or invariant at least at the critical loop levels? There are many evidences that the answer 
is yes for most of the exactly solved models but SLE is the only method so far giving us a 
machinery for rigorous proof. 

6 Conclusions 

We showed that it is possible to define a height model for a generic primary operator of 
general CFT by using the fusion matrix as an adjacency matrix. In addition it is possible to 
associate to these height models some loop models with loop weights equal to the quantum 
dimension of corresponding primary operator. 

In the critical regimes, these new loop models have some properties similar to those of loops 
corresponding to CFT. For example, the loop model corresponding to M(3, 4) has the same 
properties as the Ising model's loop representation. In the lattice level it is not completely 
clear why this should happen but it gives us a good machinery to find the loop properties of 




(5.5) 
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lattice models by studying the corresponding CFT without going to lattice level and studying 
the properties of discrete variables which most of the times are quite difficult. Investigating 
the connection with the integrable ADE models is enlightening in better understanding the 
connection with lattice models. Of course finding the CFT of the lattice model is not simple 
and so we think that our method can be useful for the more studied models with the known 
operator content. In this paper we just list some simple CFTs but as we mentioned before it 
is possible to investigate the loop models corresponding to the more general CFTs from affine 
Lie algebras to supersymmetric models and coset models. It is also interesting to generalize 
the height models to the fused cases and find the connection to the loop models with two kind 
of loops [7]. The other more interesting direction is using the method of Behrend, Pearce, 
Petkova and Zuber [3H] to find the fusion rules of boundary conformal field theory by just 
using the modular invariant partition function of the CFT and then finding the loop model 
representation by the method that we described in the paper [39J. 

Moreover we think that this method of generation of loop models can be useful to classify 
the ground state of topological quantum theories because as we already mentioned before 
the weights of the above loop models are exactly the same as quantum dimensions of the 
operators appearing in the CFT. The primary operators of CFT describe the edge states of 
topological theory and so some statistical properties of the model specially the topological 
quantum entanglement [101 HI] . 
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A Appendix: Loop Models for More General Coset Theories 

In this appendix we are going to summarize some of the possible loop models of more 
general coset theories. 

SUk su^ + S ii2) ^ : ^ ne case I = 1 was already discussed in section 3. The I = 2 is related to 
N = 1 super conformal minimal models and should have critical loops related to the loops 
of SUk(2). To get all of the critical loop representations of N = 1 superconformal models 
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one should investigate the S matrix of these models carefully. Let's see this case with more 
detail. Take k = m — 2 then the central charge of this model is c = |(1 — m ^ +2 ) ) ^ or m = 
3, 4, 5, ... and one can label the conformal weights by the triplets (j, k, I) with k = 0, 1, m, 
j = 0, 1, m — 2 and I = 0, 1, 2. We just consider the triplets with j — k + 1 being even and 
the symmetry (j, k, I) = (m — 2 — j,m — k,2 — l), so by using the S matrix given in the \27\ [28] 
we will have the following quantum dimensions 

_ sin(^)sin(^)sin(vr^) 
UAl) sin(^)sin(£)sin(f) ' ( j 

For m = 3 this theory is equal to the M(5, 4) minimal model describing three critical 
Ising model and one can find the same results as the example 1. For m = 4 one can find the 
following conformal dimensions, 



^(0,0,0) 


= 1, 




= 2, 


^(0,4,0) 


= 1, 


^(1,1,0) 


= Ve, 


d (l,3,0) 


= V6, 


d{2,0,0) 


= 1, 


^(2,2,0) 


= 2, 


^(2,4,0) 


= 1, 


^(0,3,1) 


= Vs, 


d(i,4,i) 


= 2, 


^(2,3,1) 


= V6, 


^(1,2,1) 


= 2V3. 



(A.2) 

Using equation (|3.2|) one can find the drift of CLE easily. The extension to general m shows 
that we have a primary operator (0, m — 1, 0) for odd m with the graph related to A m+ i and 
so with the following CLE drifts 

A m + 2 m + 2 

m + o m + 1 

for dilute and dense cases respectively. If we drop the constraint on the triplets then we can 

find the same equation for even m if we chose (0, m — 1,0) which is indeed possible if we 

work with modified partition function |28j . It is not difficult to show that for (m — 3,0,0) 

the corresponding graph is A m _\. These loop models can be related to the paper [6] which 

gives a lattice loop candidate for supersymmetric models. 

Finally we conjecture that for the more general case, S ^g^+K2) ^ ' we snc, uld have at 

least the following conformal curves, 

k = 4- , k = 4- , k = A- . (A.4) 

A; + 3 1 + 3 k + l + 3 y ' 

We think that these are related to critical loops of fused RSOS models introduced in [2T?] . 
The important thing to mention here is most of the loop models that we can extract with 
our method in this case are non-critical and can not be consistent. There is a method for 
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resolving this by using the coupled loop models introduced recently by Fendley [6jf which is 
based on decoupling the quantum dimension to the product of generally two different loop 
weights for fused graphs. 

SUk ~ 1 sul(N) Jl ^ ' These coset models are called Wjy algebra and describing the A n critical 
models [30], the central charge of these models is c = (N — 1)(1 — (fc + jvj(I+iv-i) )• ^ or 
these models by getting motivation from the quantum dimensions that we calculated for 
the SUk(N) one can show that for every N > 2 we have just critical loop models with the 
following CLE drift for k = 2, 

* = 4±f (A.5) 
N + 3 K ' 

It is not obvious that why the k = 2 case is an exception here but this formalism is showing 
that for k > 2 we can not find a nontrivial critical loop model. It seems that there should be 
also the same decoupling procedure as [6] however this case is totally unexplored. The above 
formula is giving the true answer for N = 3 and k = 2 which is related to three states Potts 
model. It is not difficult to see that k = 2 is related to Z(N) model and the above equation 
is in fact the same as the equation (I3.14p that we found for simplest parafermionic model. 
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